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In view of the recent LHCb measurement of AAcp, the difference between the time-integrated 
CP asymmetries in D — > K + K~ and D — > ir + ir~ decays, we perform a comparative study of the 
possible impact of New Physics degrees of freedom on the direct CP asymmetries in singly Cabibbo 
suppressed D meson decays. We systematically discuss scenarios with a minimal set of new degrees 
of freedom that have renormalizable couplings to the SM particles and that are heavy enough such 
that their effects on the D meson decays can be described by local operators. We take into account 
both constraints from low energy flavor observables, in particular D° — D° mixing, and from direct 
searches. While models that explain the large measured value for Ay4cp with chirally enhanced 
chromomagnetic penguins are least constrained, we identify a few viable models that contribute to 
the D meson decays at tree level or through loop induced QCD penguins. We emphasize that such 
models motivate direct searches at the LHC. 



I. INTRODUCTION 

Recently, the LHCb collaboration presented the first 
evidence for CP violation (CPV) in charm quark de- 
cays [1]. In particular, a difference between the time- 
integrated CP asymmetries in D — s- K + K~ and D — > 

7T + 7T~ 

A^cp, LHCb = A CP (K + K-) - A C p(7r + 

= (-0.82 ±0.21 ±0.11)% (1) 

has been reported, which is non-zero at 3.5c. This mea- 
surement is consistent at about the la level with the pre- 
vious measurement from CDF [2] , and the previous world 
average from the Heavy Flavor Averaging Group [3] . The 
new world average, combining the LHCb result with pre- 
vious measurements of Aqp(K + K~) and ^4cp(7r + 7r~) a t 
BaBar [4], Belle [5] and CDF [2], is [3] 

World Average — 

(-0.645 ± 0.180)% . (2) 

The interpretation of this measurement as a sign of 
New Physics (NP) requires a well-understood Standard 
Model (SM) calculation of this observable. Simple argu- 
ments dictate that the SM contribution to direct CPV in 
D° decays must be both CKM suppressed and loop sup- 
pressed. Concretely, the tree level decays D — > K + K~ 
and D —> ir + n~ (we implicitly include both D° and D° 
when discussing neutral D decay modes) only involve 
the first two quark generations, which cannot access the 
CP violating Kobayashi-Maskawa (KM) phase. The KM 
phase does enter into the loop-induced gluon penguin dia- 
gram for singly-Cabibbo suppressed D° decays that thus 
can provide both the required weak and strong phase dif- 
ference relative to the leading SM tree amplitude. This 
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implies that the SM prediction is loop suppressed as well 
as CKM suppressed, and the naive expectation for direct 
CPV in singly-Cabibbo suppressed D° decays is para- 
metrically given as 0{(a s /n)(V ub V* b )/{V U8 V? s )) ~ lO" 4 . 
This leads to the conclusion that the LHCb evidence of 
CPV at about the percent level is a sign of New Physics. 

A precise SM calculation, however, is difficult to ac- 
complish. Although tree level and loop level SM contri- 
butions to the quark level processes c — > uss and c — > udd 
are readily calculated, the evaluation of the hadronic ma- 
trix elements (K + K~\(uTis)(sT2c)\D°) , for example, is 
not easily performed. In the simplest approach, naive fac- 
torization, the hadronic matrix elements are "factorized" 
into (i^+|(uris|0)(i ; i:"|(sr2c)|L» ) which is formally the 
leading term in the heavy charm quark limit. As the 
charm mass is close to Aqcd, however, this approach 
suffers from large l/m c power corrections. In particu- 
lar, so-called annihilation diagrams are ignored, where 
quarks are pair-produced from the vacuum to complete 
the K or it mesons, as are long-range QCD effects such 
as final state rescattering, where constituent s quarks 
of a D — >• K + K~ decay rescatter into d quarks of a 
7r + 7r~ final state. Alternative techniques such as the 
topological diagram approach organize decay and anni- 
hilation amplitudes according to weak current insertions 
and SU (3) p light quark flavor symmetry, and then try to 
extract amplitudes and phases directly from D° branch- 
ing ratio data. 

Several recent papers have dicussed improved esti- 
mates for AA CP in the SM. In [6], a NLO QCD fac- 
torization calculation is amended by an estimate of the 
effect of certain l/rn c suppressed penguin amplitudes us- 
ing D° branching ratio data. Assuming an 0(1) strong 
phase, the authors find the SM can potentially give 
|AAcp| ~ 0.4%. As their result admittedly neglects sev- 
eral effects which could alternatively reduce or enhance 
this estimate, they conclude that the measured value of 
AAqp could be reproduced in the SM. Studies that in- 
form magnitudes and phases of D° meson decay ampli- 
tudes directly from data were performed in [7-9] follow- 
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ing a topological diagram approach. The authors of [7] 
arrive at a slightly smaller AAcp ~ —0.25% estimate, 
which remains, they highlight, more than 2a away from 
the world average. In [8, 9], the correlation between di- 
rect CPV in D — > K + K~ and h + tt~ and other D meson 
decays is emphasized as an important cross-check of the 
LHCb result. 

Even though there is large uncertainty in the SM value 
of AAcp, it is nevertheless important and exciting to 
consider the possibility that we are seeing evidence of 
NP. Literature prior to the LHCb result emphasized the 
continued fact that CPV in the charm sector is consid- 
ered an excellent probe of NP beyond the SM [10-12]. 
Among the most promising probes of CPV in the charm 
sector are observables in D° — D° mixing [13-17] and 
singly Cabibbo suppressed D decays [18-20]. In fact, 
since flavor physics observables can probe energy scales 
much higher than those directly measured, we could po- 
tentially expect that NP at the LHC would first be seen 
from its flavor effects at low energies and only later ac- 
cessed directly. 

Our goal is to investigate the possibility that NP is 
indeed responsible for the large AAqp measurement and 
to outline the corresponding NP parameter space consis- 
tent with all experimental constraints for a variety of 
NP models. Some recent work has discussed the NP 
possibility both model independently [21] and in the 
context of various concrete NP scenarios, including up- 
type flavor changing neutral currents (FCNCs), fourth- 
generation fermions, F-parity violating supersymmetry, 
and the MSSM with nonstandard sources of flavor viola- 
tion [22-26]. 

Our work differs from these previous analyses since we 
consider a much broader range of new NP possibilities 
and apply a full gamut of experimental constraints, both 
from low energy experiments and collider searches. In 
particular, we systematically discuss models with a min- 
imal set of new degrees of freedom with renormalizable 
couplings to the SM particles and are heavy enough such 
that their effects on the D meson decays can be described 
by local operators. Specifically, we consider models with 
new massive neutral gauge bosons that have flavor chang- 
ing tree level couplings to quarks, models with extended 
scalar sectors, and models where the D —> K + K~ and 
D — > 7r + 7r~ decays are modified at the loop level by 
gluon penguins. As discussed in [18, 21, 26], the loop 
induced AF = 1 chromomagnetic dipole operator (here 
and throughout, F refers to charm number except where 
noted) is, on general grounds, expected to be the least 
constrained approach for generating large nonstandard 
effects in D° meson decays. On the other hand, the 
effects of four fermion operators that are, for example, 
induced by tree level exchange of flavor changing NP 
degrees of freedom, are highly constrained by D° — D° 
mixing data. As is well known, the D° — D° constraints 
become more effective with heavier NP degrees of free- 
dom [18], leading to the expectation that almost no NP 
parameter space remains in models where four fermion 



operators are responsible for nonstandard direct CPV in 
D —> K + K~ and D — > . Our work shows to what 

extent this naive expectation holds true and identifies a 
few exceptions. We emphasize that each of the viable 
NP possibilities that can accomodate a large AAcp mo- 
tivates a further phenomenological study focusing on the 
allowed parameter space identified in this work, which we 
leave for a future study. 

In Sec. II, we review aspects of CPV in neutral D 
meson decays that are most relevant for our analysis. 
In Sec. Ill, we present the AF = 1 and AF = 2 effective 
Hamiltonians that can describe NP contributions to the 
D -> K+K- and D -> tt+tt" decays and to D° - D° 
mixing, respectively. The various NP models that con- 
tribute to CPV at tree level are discussed in Sec. IV, 
while the NP models that contribute at loop level are 
discussed in Sec. V. We conclude in Sec. VI. Technical 
details about hadronic matrix elements and renormaliza- 
tion group running, as well as a collection of loop func- 
tions can be found in the appendices. 



II. CP ASYMMETRIES IN NEUTRAL D 
MESON DECAYS 

The neutral D meson mass eigenstates D\ and D 2 are 
linear combinations of the strong interaction eigenstates, 
D° and D° 



\D li2 ) =p\D°)±q\D°) 
The factors q and p are given by 
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where M 12 and Ti2 are the dispersive and absorptive part 
of the D meson mixing amplitude. CP violation in D 
meson mixing is signaled by \q/p\ ^ 1 or cf> = Arg(g/p) ^ 
0. 

The normalized mass and width differences, x and y, 
in the neutral D meson system are given by 
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(6) 



where the lifetime of the D° mesons t = 1/T = 0.41 
ps [10]. _ 

The time integrated CP asymmetry in the decay of 
neutral D mesons to a final CP eigenstate / = K + K~, 
7r + 7r _ is defined as 



Acp(f) 



T(D° -» /) - r(F° -> f) 

A m + A i + A d ( 7 ) 
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The time integrated CP asymmetry receives contribu- 
tions from CPV in mixing A m , CPV in interference of 
decays with and without mixing A 1 , and from CPV in 
the decay itself A'j. The "indirect" CP asymmetries A m 
and A % are approximately independent of the final state 
and depend only on D° — D° mixing parameters 
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where 7/q P is the CP parity of the final state. The 
"direct" CP asymmetry A'j is instead sensitive to the 
final state. The decay amplitudes of singly Cabibbo 
suppressed D meson decays A(D° — »• /) = At and 
A(D° ->• /) = Af can be written as [18] 



A f = A\ 1 



A, 



^(l + r/e^/-*/) 



J 



(10) 
(11) 



where AJ is the dominant singly Cabibbo suppressed tree 
level amplitude, which can be taken real by convention, 
and rf is the relative size of subleading ("penguin") am- 
plitudes. With respect to the tree amplitude, the pen- 
guin amplitudes can have a relative weak phase <j)f and 
a relative strong phase 8f. 

Under the assumption that rf is small, one arrives at 
the following expression for the direct CP asymmetry 



A' 1 



f 



2r/ sin 5f sin <f>f 



(12) 



The difference between the time-integrated CP asymme- 
tries in D — > K + K~ and D — » ir + ir~ measured by LHCb 
is given by [1] 



AA 



CP 



.4 



K+K- 



A d v+ ^ + ^±(A m + A l ) , (13) 



where A(t)/r = (9.8 ± 0.9)% is a small difference in the 
average decay times of the D° mesons in the K + K~ and 
7r + 7r~ sample [1]. Given the existing bounds on the in- 
direct CP asymmetries [3], the LHCb measurement of 
AAcp is an excellent approximation of the difference in 
the direct CP asymmetries. 

As already mentioned in the Introduction, charm CPV 
in the SM is strongly Cabibbo suppressed. Furthermore, 
in the SM, direct CP violation in D -> K + K~ and D 
7r + 7r~ decays comes from the interference of the tree level 
contribution with a loop suppressed penguin amplitude 
and correspondingly, r f - 0(a s /ir){V ub V* b )/(V us V* s ) ~ 
10~ 4 . Even though the weak phase of the SM penguin is 
large (7 ~ 70°) and assuming a maximal strong phase, 
a naive SM estimate for AAqp is therefore smaller than 
the global average by at least an order of magnitude. 

Sizable direct CP asymmetries in the D — !• K + K~ 
and D —> 7r + 7r~ decays are only possible in the SM if 



the relevant hadronic matrix elements are strongly en- 
hanced [27]. Despite several recent studies [6-8], it re- 
mains unclear to what extent such an enhancement is 
present and whether the value of A^4cp measured by 
LHCb can be explained within the SM. 

In the following we investigate the possibility that the 
measured AAcp is due to New Physics. 



III. EFFECTIVE HAMILTONIAN APPROACH 

A. AF = 1 Effective Hamiltonian 

In the New Physics frameworks discussed below, con- 
tributions to the singly Cabibbo suppressed D — > K + K~ 
and D — > tt + tt~ decays can be described by the following 
effective Hamiltonian 



P 2 — 1 



J2 [ c i °i °i 



h.c. 



(14) 



where X p — V cp V* , and the operators 0\ ' are given by 
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(15c) 
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(15i) 
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,m c ua^(l + 75 )c^G- u , (15k) 



(uP L s)(sP L c) , 
{u a PLSp)(spP L c a ) , 
{ua^P L s){sa^P L c) , 



O k t I = {u a (T ll . v P L sp){spa tlv P L c a ) 



(151) 
(15m) 
(15n) 
(15o) 
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The index q runs over all active quark flavors, the index 
p runs over all active down type quark flavors, a and 
P are color indices (that are implicitly summed over), 
e q is the electric charge of the quark q, (V ± A) refers 
to the Dirac structures 7 M (1 ± 75), Pr,l = 5(1 ± 75) 

and = \idn1v ~ Ivlp)- The operators 0,[ 1){p} are 

obtained from O^ 1 '^ by replacing 75 — > —75. 

The operators Oj 1 ^ are the so-called current-current 
operators. In the SM, tree level W exchange generates 
at the matching scale the Wilson coefficient C-[ ~ 
Gf/\/2. The QCD penguin operators 5 6 and the 

chromomagnetic operator O^J are first generated at 
O(ot s ) and proportional to V u bV* b . The chromomagnetic 
operator is proportional to the charm quark mass but can 
be chirally enhanced by v/m c from NP. The QED pen- 
guin operators ofy 9 10 are also proportional to V u bV* b . 
They are of 0(a) and negligible in the D — > K + K~ 
and D — > ir + ir~ decays in the SM. In the NP models 
discussed below that have tree level contributions to the 
D —> K + K~ and D — > it + it~ decay amplitudes, however, 
the QED penguin operators can be relevant. The scalar 
operators Ogl S2 become important in the context of the 



2HDM discussed in Subsec. IV E and the scalar octet dis- 
cussed in Subsec. IV F. The tensor operators 0^\ T2 do 
not contribute to D — > K + K~ and D — » 7r + 7r~ decays 
in naive factorization. We consider them nonetheless, 
because they mix with the scalar operators under renor- 
malization group running. 

The ratio 77 that enters the expression for the direct 
CP asymmetry Eq. (12) can be written as a function of 
the Wilson coefficients appearing in Eq. (14). We use 
the results from [18] for the hadronic matrix elements 
that are obtained using naive factorization for 6 and 

QCD factorization [28, 29] for O^J . The matrix elements 
obtained in naive factorization are formally the leading 
terms in an expansion in a s and Aqcd/^c [28, 29]. In 
the case of D meson decays, however, and as mentioned 
in the Introduction, it is known that power corrections, in 
particular annihilation contributions, which are formally 
suppressed by l/m c , can be equally important [6, 8, 18]. 
The naive factorization results can therefore only be con- 
sidered as rough estimates and in our numerical analysis, 
we will allow for enhancements up to a plausible factor 
of 3 [6, 18]. For our analysis, we extend the results for 
the hadronic matrix elements given in [18] by including 
the QED penguin and scalar operators (see Appendix A 
for details). We find 
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where iV c = 3 is the number of colors and p = s, f — 
K + K~ for the D — > K + K~ decay and p = d, f = tt + it~ 
for the D — > ir + ir~ decay. The chiral factors Xf are 
approximately given by 

^ 2m K v , 17 n 

Xk+k- — j Xit+tt- — / ; 7 ! (J- 'J 

m c m s m c (m d + m u ) 

with all quark masses evaluated at the scale of the D 
meson n ~ mo ~ 1.8 GeV. All the Wilson coefficients 
in Eq. (16) are evaluated at this scale. We use LO renor- 
malization group running for 10 and Co\ S2 ti T2 as 

well as for C^J to evolve the Wilson coefficients from 
the high matching scale, where NP degrees of freedom 
are integrated out, down to fi ~ rri£>. The correspond- 
ing anomalous dimensions are collected in Appendix B. 
We do not include 2-loop mixing between and the 
other Wilson coefficients. In view of the large uncertain- 
ties in the evaluation of the hadronic matrix elements, we 
consider this approximate approach to be fully justified. 
While there are no strong phase differences between 



the several operator matrix elements in the naive factor- 
ization approach, they can be generated by large power 
corrections or final state interactions. Throughout this 
work, we will assume 0(1) strong phase differences, fol- 
lowing [6, 18]. 

B. AF = 2 Effective Hamiltonian 

In the models discussed below, the most important fla- 
vor constraints come often from D° — D° and K — K 
mixing. New Physics contributions to meson mixing can 
be described by the effective Hamiltonian 

In the case of D° — D° mixing, the most important op- 
erators for our analysis are given by 

Of )D = {u al ^P L c a ){ui3Y P L Cfj) , 
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6f )D = {u a ^^P R c a ){u^P R ci3) , 

6f ] ° = (Uai 5 flC«)(%P R ^) . (19) 

In the case of K — K mixing, the operators most relevant 
for our analysis are 

of )K = {d a 7 M P L 8 a )(d f n»P L 8f ) ) , 

Of )K = {d a 7 M P R 8 a )(df,<fP R 8f l ) , 

0f )K = (d a P L s a )(dpP RSp ) , 

Of )K = (d a P L s p )(dpP R s a ) . (20) 

In the above expressions, P r .l = j(l± 75) and a, f3 are 
color indices (that are implicitly summed over). 

The Wilson coefficients C\ are again obtained by in- 
tegrating out the NP degrees of freedom at a scale of the 
order of the mass of the new particles. Using renormal- 
ization group evolution [30, 31], these coefficients are sub- 
sequently run down to the low scale where the hadronic 
matrix elements [32-34] are given. Combining Wilson co- 
efficients with the hadronic matrix elements gives the NP 
contribution to the dispersive part of the mixing ampli- 
tude M 12 . 1 In the case of Z?° — D° mixing, the SM contri- 
butions to neither the dispersive part nor the absorptive 
part of the mixing amplitude can be predicted reliably as 
they are dominated by long distance effects [35, 36]. In 
our numerical analysis, we allow the long distance contri- 
butions to vary in the ranges M^ 2 D £ [-0.02,0.02] ps^ 1 
and r^ 2 D G [-0.04,0.04] ps" 1 [37], such that by them- 
selves they can saturate the experimental values. We 
apply the most recent averages and la errors of the ex- 
perimental results on the D° — D° mixing parameters [3] 

x = (0.63ig;g)% , y = (0.75 ± 0.12)% , 

1^1=0.891°;^, 0= (-10.lt!;*) , (21) 

at the 2a level throughout our analysis. 

We note that in many of our minimal field content 
scenarios, the NP vertices used in the AF = 1 operators 
are also used for the AF = 2 operators, leading to a phase 
relation 2<f>p = x — 4>f=i between the CPV for D° decays 
and the CPV for D° — D° mixing. This relation implies 
that the D° — D° mixing constraint is best satisfied by 
eliminating the CPV in D° - D° mixing and saturating 
the D° — D° mixing transition amplitude. On the other 
hand, in non-minimal constructions, this phase relation 
could be different, possibly making the null observation 
of CPV in D°—D° mixing the more restrictive constraint. 

IV. NEW PHYSICS CONTRIBUTIONS AT 
TREE LEVEL 

We concentrate on New Phyiscs models where the new 
degrees of freedom are heavy enough such that their ef- 
fects in low energy observables can be reliably described 



1 The absorptive part of the mixing amplitude Ti2 is not sensitive 
to new short distance dynamics. 
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u 

(a) (b) 

FIG. 1. Example tree level Feynman diagrams that con- 
tribute (a) to the D — > K + K~ and D — > 7r + 7r~ decay am- 
plitudes and (b) to D° - D° mixing in the cases of a flavor 
changing Z, Z' and heavy gluon G' . 

by the local operators introduced in Sec. III. We do not 
consider scenarios with very light mediators, which is be- 
yond the scope of this work. Moreover, we focus on mod- 
els where the new degrees of freedom have renormalizable 
couplings to SM degrees of freedom. 

In this section, we analyze models where at most one 
new field is added to the SM. We first discuss exten- 
sions of the SM in which a massive neutral gauge bo- 
son leads to tree level contributions to the D — > K + K~ 
and D — > 7r + 7r _ decay ampitudes. We consider a flavor 
changing coupling of the SM Z boson in Subsec. IV A, 
a flavor changing Z' in Subsec. IV B and a flavor chang- 
ing heavy gluon in Subsec. IV C. We also comment on 
the possible effects of a new charged gauge boson in Sub- 
sec. IV D. Then, we analyze models with extended scalar 
sectors, namely a 2 Higgs doublet model with Minimal 
Flavor Violation in Subsec. IV E, a model with a scalar 
octet in Subsec. IV F and a model with a scalar diquark 
in Subsec. IV G. 

Models that contain more than one non-SM particle 
and where NP contributions to the D meson decays are 
first generated at the one loop level are discussed in 
Sec. V. 



A. Flavor Changing Z 

We consider a flavor changing coupling of the SM Z 
boson to the right-handed charm and up quark 

Ant = X cu c R ^u R Z^ + h.c. , (22) 

where X cu is a complex parameter. A complementary 
setup, where flavor changing couplings involving the top 
quark generate an effective c — > u transition at the loop 
level by a double flavor flip c — > t — > u, is discussed in [26]. 
Flavor changing Z couplings can appear in various sce- 
narios [38], for example in models with non-sequential 
generations of quarks [39] and also in models with ex- 
tra U(l) gauge symmetries [40], or can be loop induced 
as in SUSY models. In the absence of SU(2) L breaking 
sources, the c R u R Z coupling has the form of a charge ra- 
dius interaction and vanishes for q 2 — > 0, where q is the 
momentum of the Z boson. The dominant contribution 
to the coupling X cu is therefore in general expected to 
be proportional to t) 2 /A^ P , where v is the Higgs vacuum 
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of 3. The red (dash-dotted) region is excluded by the con- 
straints from D° — D° mixing. The D° — D° constraints 
are minimized for Arg(X cu ) — 0, tt/2, tt, 37r/2, where 
constraints from CPV in D° — D° are not effective and 
the dominant constraint comes from the normalized mass 
difference x. Indeed, if no enhancement of the hadronic 
matrix elements is assumed, sizeable NP effects in AAqp 
are only compatible with D a — D° mixing in a small cor- 
ner of parameter space with Arg(X cu ) ~ 7r/2, 37t/2. Still, 
barring the finetuned situations Arg(X c „) = 7r/2,37r/2, 
sizeable NP effects in A^4cp a l so imply indirect CPV 
in D° — D° mixing close to the current experimental 
bounds. The required size of the flavor changing cou- 
pling v 2 /Afj P oc \X CU \ ~ 10~ 4 points towards a NP scale 
of Anp < few x 10 TeV, where this coupling is generated 
and not necessarily within the immediate reach of direct 
searches. 



B. Flavor Changing Z' 



FIG. 2. Regions in the |X CU | - Arg(X cu ) plane compatible 
with the data on AAcp, wa at the la level in the model with a 
flavor changing Z. The green (solid) band corresponds to the 
expressions for the decay amplitude in naive factorization, the 
blue (dashed) band assumes an enhancement of the hadronic 
matrix elements by a factor of 3. The red (dash-dotted) region 
is excluded by the D° — D° mixing constraints. 



expectation value (vev) and Anp is the NP scale where 
the flavor changing Z coupling is generated. 

As shown in diagram (a) of Fig. 1, the X cu coupling 
leads to tree level contributions to the Wilson coefficients 



Cf\ &!p and C, 



(i) 



a 



(i) _ 



1 9 X* cu 
3 2c w 4M| 



(i) _ 



cu 



a 



(i) 



Z 9CW AMl 

2 gs 2 w x* 



3 c w 4M| 



(23) 



The flavor changing cuZ coupling also inevitably gener- 
ates tree level contributions to D° — D° 



mixing 



2)D 



2M| 



(24) 



If the Z boson has flavor changing couplings to left- 
handed quarks, 1-loop contributions to e'/e would also 
be generated. In order to avoid this constraint, we re- 
strict ourselves to the criirZ coupling. 

In Fig. 2, we show the regions in the \X CU \ - Arg(X cu ) 
plane that are compatible with the range for A^4cp in 
Eq. (2) at the lcr level. The green (solid) band is ob- 
tained using the expressions for the decay amplitude in 
naive factorization. The blue (dashed) band assumes an 
enhancement of the hadronic matrix elements by a factor 



Next, we consider a leptophobic massive Z' gauge bo- 
son with tree level flavor changing couplings to right- 
handed up and charm quarks. Models of this type can 
be easily constructed if the Z' couples with SM degrees of 
freedom through higher dimensional operators [41]. Fla- 
vor changing couplings of a Z' can also arise, for example, 
in models with family non-universal couplings [40, 42]. 
We parameterize the interactions of the Z' in the follow- 
ing way 



C 



inn = gL^Tu^Z'^ + g u u R ^u R , 
+g L dl^dlZl + g d d Rl »d R Z^ 
+X c R ^u R Z' + h.c. , 



(25) 



where the flavor universal couplings g u , g& and gi, are free 
real parameters and the (small) flavor changing coupling 
X cu is a free, complex parameter. We restrict ourselves to 
a tree level c — > u coupling and do not consider t — > c and 
t — > u couplings that could induce the c — > u transition 
at the loop level. 

Depending on whether the flavor conserving couplings 
of the Z' are to left-handed or right-handed quarks, the 
flavor changing cuZ' coupling can induce tree level con- 



tributions to the Wilson coefficients and C, 



(i) 





{g u + 2gd) 

3 


X* 

CU 

4M§, 


c« 


2 (5« - gd) 
3 


X* 

cu 

AM\, 




9lX* u 
~ AMI, 





If we assume g. 



(26) 



i) 



and C L 
the D 



(i) 



gd, then only contributions to C 3 
are generated. The contributions from C^l 9 to 
K + K~ and D — > ti + it~ decay amplitudes are 
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FIG. 3. The M z > — \X CU \ plane, setting Arg(X cu ) = n/2. In the left plot, AAcp is evaluated in naive factorization and in the 
right plot, we allow for an enhancement by a factor of 3. Along the solid green and blue lines the NP contributions to AAcp 
match the world average. The different green and blue lines correspond to different common choices of the flavor conserving 
couplings g u ~ gd ~ 9l as indicated. The black dashed line shows the constraint on the flavor conserving coupling from dijet 
searches. The red (dash-dotted) region is excluded by D° — D° mixing constraints. 



color suppressed. The contribution from Cg is helicity 
enhanced. 

The flavor changing cuZ' coupling also induces tree 
level contributions to D° — D° mixing 



(2)D 



(x* c 



2Mf, 



(27) 



We highlight the following point: as both the NP con- 
tributions to the AF — 2 mixing amplitude and the 
AF = 1 decay amplitudes are described by dimension 
6 operators, they decouple with the NP mass squared. 
Yet while the AF = 2 amplitude is obviously propor- 
tional to the square of the flavor changing coupling, the 
AF = 1 amplitude is linearly proportional in this cou- 
pling. Correspondingly, the constraint from D° — D° 
mixing becomes more effective with heavier NP mass. 
Analogous arguments hold in all the other NP scenarios 
discussed in this work. 

If the Z' boson has flavor changing couplings to left- 
handed quarks, tree level contributions to K — K mixing 
would be unavoidably generated since the clu^Z' and 
slcIlZ' couplings are related by the CKM matrix due to 
SU{2)l invariance. Since constraints coming from K — 
K mixing are considerably stronger than those coming 
from D° — D° mixing, we restrict ourselves to the crUrZ 
coupling. As a result, 1-loop contributions to e'/e that 
can lead to constraints are also absent. 

The dijet searches at hadron colliders set additional 
constraints on the model. In this paper, we consider the 
searches at UA2 [43], CDF [44] and CMS [45]. The UA2 
collaboration probed the light dijet mass region, from 



130 GeV to 300 GeV, while the CDF search covers a 
dijet mass range from 260 GeV to 1.4 TeV. A higher 
mass range, 1.0 — 4.1 TeV, is probed by the CMS exper- 
iment. There are other dijet searches from the D0 [46] 
and ATLAS experiments [47]. The D0 collaboration an- 
alyzed 109 pb _1 of data, however, while CDF analyzed 
1.13 fb~ , indicating the D0 bound is less competitive 
than the one from CDF. The ATLAS bound is expected 
to be comparable with the CMS bound since both of the 
experiments analyzed 1 fb _ of data. 

We simulate the Z' production using MadGraph 5 [48] , 
and the width of the Z' is calculated with CompHEP [49], 
varying g u = gd = 9l- We compare the simulated cross 
section with the limit on dijet production from UA2 
(Fig. 2 of [43]), CDF (Table I of [44]) and CMS (Ta- 
ble 1 of [45]). In calculating the bound, we ignore X cu 
since it is at least one order of magnitude smaller than 
9u, gd and g L . 

The plots in Fig. 3 show the M z > — \X CU \ plane, setting 
Arg(X cu ) = 7r/2. In the left plot, AAqp is evaluated in 
the naive factorization approach, while in the right plot 
we allow for an enhancement of the hadronic matrix ele- 
ments by a factor of 3. Along the green (solid, left plot) 
and blue (solid, right plot) lines, the NP contributions 
to AAqp match the world average. The different green 
or blue lines correspond to different choices of the flavor 
conserving couplings g u = 9d = 9l = 0.1, 0.2, or 0.5. The 
region below the black (dashed) line requires a flavor con- 
serving coupling that is excluded by dijet searches. The 
red (dash-dotted) region is excluded by the constraints 
from D° — D° mixing. The choice Arg(V c „) = ir/2 cor- 
responds to a maximal phase for the NP contributions 



to the D — > K + K~ and D —> ir + ir~ decays while si- 
multaneously minimizing the constraint from D° — D° 
mixing. Choosing Arg(X cu ) = n/3 would lead to C(l) 
phases both in the decays and in D° — D° mixing, and 
the corresponding constraint would be more stringent by 
a factor of ~ 2. As expected, the D° — D° constraint 
becomes more effective with larger Z' mass. 

We observe that even allowing for an enhancement in 
AAqp by a factor of three, the D° — D° mixing constraint 
in combination with dijet searches rules out a Z' as a 
tree level NP explanation for the measured AAqp- We 
do not consider Z 1 masses below 100 GeV, which would 
be constrained from Z — Z' mixing [50], but the exact 
constraints would be model dependent and are beyond 
the scope of this work. 



C. Flavor Changing Heavy Gluon 

Heavy color octet vector bosons that have flavor chang- 
ing couplings to SM quarks can arise, for example, 
in models with warped extra dimensions [51] and also 
in models of axigluons with family non-universal cou- 
plings [52]. We concentrate again on a direct tree level 
coupling between right-handed up and charm quarks 

Ant = g L uirT a u l L (G% + g u u^T a u l R {G') a 



butions to D° — D° mixing 



+g L d l L rT a dl(G% + g d d Rl »T a d R {G') 
+X cu c R ^T a u R {G')l + h.c. . 



/ \ a 



(28) 



The flavor universal couplings g u , g& and g^ are free real 
parameters: the (small) flavor violating coupling X cu is 
a free, complex parameter. 

The heavy gluon can generate tree level contributions 



to ^3,4,5,6 an d C'lio 



_ Cgg + 2ff rf ) X* 



C 



8M G , 



Mi) _ ~ l Mi) 

° 3 ~vT 4 ' 



Mi) _ %9u ~ 9d) Mi) Mi) _ 

9u + £g<i -<v c 



Mi) "/--I 

8M* ' °5 - AT °6 



Nr. 



(29) 



To generate and Cg 1 ^, the RH flavor conserving cou- 
pling is required, while for Cgg the LH flavor conserving 
coupling is required. If we assume g u = gd, Cgl are ab- 
sent. The Wilson coefficients C35 9 are color suppressed, 
and furthermore, their contributions to the decay ampli- 
tudes are color suppressed. The contributions to the de- 
cay amplitudes from and Cg 1 ^ are helicity enhanced. 
Due to the presence of Cg , we expect slightly larger NP 
contributions to AAqp than the Z' scenario. 

Similar to the models discussed before, the flavor 
changing cuG' coupling also leads to tree level contri- 



(2)D _1-N C (X*J 2 



c[ 2)D = 



2N n 



2Ml, 



(30) 



Compared to the Z 1 case, the contribution to D° — D° 
mixing is suppressed by a factor |(1 — N C )/(2N C )\ = 1/3. 
As discussed for the Z' case, we do not consider flavor 
changing couplings to left-handed quarks to avoid the 
stringent constraints from K — K mixing and e'/ e - 

Collider constraints come again from dijet searches at 
the hadron colliders, and we evaluate them in a similar 
fashion to the Z' case. Additional constraints arise from 
recent results on four jet searches at the LHC, searching 
for pair production of dijet resonances. Pair production 
of the heavy gluon is fixed by QCD to good approxi- 
mation and the production cross section depends only 
on the G' mass [53]. The mass range from 100 GeV 
to 200 GeV is covered by an ATLAS search using 34 
pb^ 1 [54] while a CMS search using 2.2 fb" 1 [55] starts 
at 320 GeV. The intermediate region from 200 GeV to 
320 GeV is not considered in the CMS search because 
of the multijet trigger turn-on effects on the QCD back- 
ground fit curve, which well models the QCD multijet 
background above 320 GeV. We simulate pair produc- 
tion of the G' resonance using Madgraph 5. We find that 
a G' mass in the ranges 100 GeV < Mqi < 200 GeV is ex- 
cluded. Assuming the signal acceptance to be 3% (13%), 
which corresponds to the lowest (highest) signal efficiency 
found in [55] , we find that a G' mass from 320 GeV up to 
720 GeV (1000 GeV) is excluded. We conservatively use 
the stronger bound derived from an acceptance of 13% 
in the following discussion. 

The plots in Fig. 4 show the Mq' — \X CU \ plane with 
Arg(A ctl ) = 7r/2 in order to minimize the constraint 
from D° — D° mixing. In the left plot, A^4cp is eval- 
uated in the naive factorization approach, while in the 
right plot we allow for an enhancement by a factor of 
3. The NP contributions to AAqp match the world av- 
erage along the green (solid, left plot) and blue (solid, 
right plot) lines. The different green or blue lines cor- 
respond to different choices of the flavor conserving cou- 
plings g u = 9i = gh = 0.1, 0.2, or 0.5. The region below 
the black (dashed) line requires a flavor conserving cou- 
pling that is excluded by dijet resonance searches. The 
red (dash-dotted) region is excluded by the constraints 
from D° — D° mixing. The vertical bands are excluded 
by the dijet pair searches at ATLAS (left band) and CMS 
(right band). 

Because of different 0(1) factors in the D — > K + K~ 
and D — > 7r + 7r~ decay amplitudes and D° — D° mixing 
compared to the Z' , a heavy gluon appears slightly better 
suited to generate nonstandard effects in AAqp ■ Yet only 
after allowing for an enhancement in AAqp by a factor 
of 3 can the combined constraints from D° — D° mixing 
and dijet searches be made compatible with the measured 
AAcp. The corresponding corner of parameter space is 
characterized by light G' masses Mqi < 300 GeV. As 
the currently available results for dijet pair searches do 
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FIG. 4. As Fig. 3 but in the model with flavor changing heavy gluon. The shaded vertical bands are excluded by dijet pair 
searches at ATLAS (left band) and CMS (right band). 



not exclude the range between 200 GeV and 320 GeV, a the SM couplings gV c i,V* b in order to generate a sizable 

heavy gluon cannot be ruled out as a possible NP expla- A^cp. Such a W would then lead to unacceptably large 

nation of the observed AAcp. If a heavy gluon is indeed NP contributions to B —> D and B — > it transitions, and 

responsible for the large value of AAqp, indirect CPV in thus we will not consider this scenario any further. 
D° — D° mixing is also expected to be close to the current 
experimental bounds. 

E. Two Higgs Doublet Model 



D. Charged Vector Boson 

We consider a new vector boson with charge ±1 that 
couples to right-handed up and down type quarks 

A„t = g R Vif d l B u J R W'- + h.c. , (31) 

where is a unitary mixing matrix, the analog of the 
CKM matrix in the right-handed sector. One possibil- 
ity to introduce such a W gauge boson is through an 
additional SU{2)r gauge group [56]. Yet as long as the 
coupling structure in Eq. (31) is realized, the exact imple- 
mentation of the W is of no relevance for the following 
discussion. 

Tree level exchange of the W gives contributions to 
the current-current Wilson coefficient C^ p . The corre- 
sponding operator 0^ p has the same hadronic matrix 

element as the SM operator oj 1 ^. Therefore, tree level 
exchange of the W cannot generate a direct CP asym- 
metry because there is no strong phase difference with 
respect to the LO SM contribution. 

Loop level contributions, i.e. gluon penguins with W' 
loops, have a structure that is analogous to the SM pen- 
guin contribution. If we assume that the mass of the W' 
is larger than the mass of the SM W boson, then the 
couplings gnV^V^* have to be considerably larger than 



One of the simplest extensions of the SM scalar sector 
is the 2 Higgs doublet model (2HDM) (see [57] for a re- 
cent review). The most general couplings of the 2 Higgs 
bosons to SM fermions read 

Ant = Y U QUH U + Y d QDH d + Y t LEH d (32) 

+ X U QUH\ + X d QDHl + XiLEHl + h.c. , 

where we suppress the flavor indices on the quark fields 
and the coupling matrices X and Y. Generically, the 
neutral components of both Higgs bosons acquire a vev 
(H® d ) = v Ut d and the fermion mass matrices receive con- 
tributions from both Y and X. As Y and X are indepen- 
dent matrices in flavor space, the couplings in Eq. (32) 
lead to flavor changing neutral Higgs vertices at tree level 
and are strongly constrained from meson mixing observ- 
ables. 

The most effective way to control flavor changing ef- 
fects in the 2HDM is the Minimal Flavor Violation 
(MFV) ansatz [58, 59]. Applied to the 2HDM, the MFV 
assumption states that the "wrong" Higgs couplings X 
can be expanded in powers of the Yukawa couplings Y 

X u = £uY u ~\~ £ U Y U Y^Y U -\- £ u Y d Y^Y u -f- . . . , 

X d = e d Y d + e' d Y u YjY d + e'^Y d YjY d + ... , (33) 

where the are free complex parameters. For simplicity, 
we will also assume Xg — tiYe. 
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FIG. 5. Example Feynman diagrams that contribute to (a) the D — > and D — > 7r + 7r decay amplitudes, (b) -D — D° 

mixing, (c) K — K mixing, (d) the B + — ¥ tv decay and (e) the Bd — > X a y decay in the discussed 2HDM with MFV. 



We stress that in this particular 2HDM, there exists no 
preferred basis for the 2 Higgs doublets. The Ci param- 
eters as well as tan/3 = v u /vd are basis dependent and 
none of them separately are actually physical parame- 
ters (see [60] for a detailed discussion). In the following 
analysis, we fix a basis by setting ed = 0. In this basis, 
a large tan /3 can be approximately identified with the 
basis invariant enhancement of the coupling of the right 
handed strange quark to the charged Higgs with respect 
to its SM Yukawa coupling. 

We now investigate the parameter space of this 2HDM 
with MFV in a basis with ed = 0. We work in the regime 
of large tan/3 and assume e", e' d 1 in order to ensure 
that tree level FCNCs are under control. Furthermore, 
we allow the parameters e u ,e' u and e d as well as et to 
be 0(1). As we will see below, in this region of pa- 
rameter space, sizable nonstandard contributions to the 
D — > K + K~ amplitude can arise, while the most impor- 
tant constraints can be kept under control. 

In this scenario, tree level charged Higgs exchange, as 
shown in diagram (a) of Fig. 5, gives the dominant NP 
contribution to the D — > K + K~ decay. For large values 
of tan /3 and assuming e u ~ O(l) we find 2 



C 



(i) 

51 



m c m s 



tan/ 



V V* 



1 



..tan/9 Ml 



(34) 



with v 2 = v 2 . +Vd = 174 2 GcV 2 , and e s = e^'y 2 . For large 
tan/3 and e" d ~ 0(1), we find e s ~ 10 -3 , such that e s 
only becomes relevant for extremely large tan /3. The pa- 
rameter e„ simultaneously lifts the 1 / tan /3 suppression 
of the crslH + vertex and provides a source of CPV. 



2 Note that this expression as well as the ones given below are not 
basis invariant. They only hold in bases where the same condi- 
tions on the ti and tan fj hold, under which they were derived. 



The 2HDM with MFV have been thoroughly studied 
in the literature (see [59, 61-64]) and various constraints 
have been identified. In the following, we discuss the 
most important constraints: 

(i) Direct searches at LEP for a charged Higgs give the 
bound M H ± > 80 GeV [65]. 

(ii) The B + — > t + v and K + — > decays are 
known to be important low energy probes of extended 
Higgs sectors. Combining the experimental results from 
Belle [66, 67] and BaBar [68, 69] on Br{B tv) with 
a conservative SM prediction based on \V u b\ = (3.89 ± 
0.44) x 10" 3 [10] and f B + = 196.9 ± 8.9 MeV [70] 

Br{B -» tj/) sm = (1.04 ± 0.25) x 10" 4 , (35) 

we find 

Br(B -> tv) sm 

In our 2HDM, tree level exchange of a charged Higgs (see 
diagram (d) in Fig. 5) results in 



Rbtv 



1 



tan/3 



Mg± 1 + e& tan /3 ti 



(37) 



where we defined if, = (" d vi- For large tan/3 and (" d ~ 
0(1), one typically has q, ~ 10 -2 — 10 _1 . The factors ib 
and ei keep the bRU^H^ 1 and tbMi j H ± couplings small 
for large tan /3 and the experimental constraint from 
B + —> t + v is easily avoided, provided sign(e£) = — 1. 
In fact, for such choices of parameters, the ~ 2tr discrep- 
ancy between the SM prediction and the experimental 
result for Br(_B — > tv) is reduced. 

For the observable R123 [71], which is sensitive to 
charged Higgs contributions to the K —} (j,v decay, we 
find 



tan/3 



1 



e s tan /3 ei 



(38) 
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From the experimental side, one has [71] 
R123 = 0.999 ± 0.007 , 



(39) 



which, in our framework, leads only to constraints for 
extremely large tan /?. 

(iii) Charged Higgs loops lead to contributions to 
D° — D° and K — K mixing (see diagrams (b) and (c) 
of Fig. 5). In the considered scenario, we find the follow- 
ing dominant NP contributions 



(2)D 



-1 < (V CS V: S ) 2 tan 4 /? 



* \2 



: tan/?| 4 



(40) 



(a)J r -1 m\ (V ts V t * d ) 



4M 2 H± 



(jit^hiixt) + \i t \ 2 h 2 (x t ,x w ^ , (41) 



with x t = m 2 /M 2 H± , x w = M^/M^ ± , e t = e u + e' u yf, 
and the analytical expressions for the loop functions hi 
and h 2 are given in Appendix C. The expression for C^ D 
is suppressed by four powers of the strange quark mass 
and is only relevant for very large tan/3. 3 The Wilson 
coefficient C^ K arises from box diagrams including both 
one W+ and one H insertion as well as diagrams with 
two H+ insertions. Its contribution to kaon mixing can 
be relevant for i t — 0(1). Yet even for |et| ~ 0.5, the 
constraint on the charged Higgs mass is as low as the 
bound from direct searches M H ± > 80 GeV. 

(iv) As shown in diagram (e) of Fig. 5, the Bd — > X s j 
decay also receives 1-loop charged Higgs contributions. 
As is well known, the good agreement of the experimen- 
tal data and the SM prediction of its branching ratio 
leads to the constraint M H ± > 300 GeV in a 2HDM of 
type II [72]. This bound, however, does not apply in the 
model considered here. For the ratio of the b —> am- 
plitudes in our 2HDM with MFV model to the 2HDM of 
type II, we find, to a good approximation, 



A(6^s7)mfv e t tan/3 2 
Mb -> S7 )„ = l + e> tan/3 + N ^ ' 



(42) 



where y(xt) is a function of xt = m 2 /M H ± and is 0(1). 
We point out that because of the it and it, factors, the 
b — > «7 amplitude can be complex, which relaxes the 
Bd — s- X s j constraint considerably [73]. Still, charged 
Higgs masses as low as the constraint from direct searches 
require i t — e u + e' u y 2 <C 1, which implies a considerable 
amount of fine tuning. 

(v) Complex parameters also lead to contributions 
to electric dipole moments (EDMs) (see [61, 62]). While 
a detailed study of EDMs is beyond the scope of this 
work, we mention that in the studied framework, the 
contributions to the EDMs depend on the parameters it 



3 The corresponding contribution from a bottom quark loop is 
strongly Cabibbo suppressed and turns out to be much smaller. 



and it,- Therefore, by allowing for a certain amount of 
fine tuning, EDM constraints do not exclude sizable CP 
violating effects in the D — > K+K~ decay that depend 
mainly on e u . 

(vi) Charged Higgses are also constrained by bounds 
on the branching ratio of the top quark decay t — > H + b. 
ATLAS and CMS obtained bounds at the level of Br(< 
H+b) < 5% under the assumption Br(H~ — » tv) — 
1 [74]. CDF and D0 also give bounds considering the 
H —> cs final state. These bounds are at the level of 
Br(t -> H+b) < 10% - 20% [75, 76]. In our setup, the 
ratio of the t — > H+b and t — > Wb branching ratios is 
given by 



?tHb _ 1 - {m H /m t ) 2 

Ftwb 1 - (m w /m t ) 2 m^m 2 
ml tan 2 /3 



m z t \itY 



|1 + e 6 tan/3| s 



(m 2 - m 2 H ) 



+4Rc 



it tan (3 
1 + if, tan j5 



(43) 



For the parameter choices detailed below, we find the 
charged Higgs branching ratio into tv does not exceed 
10% due to the strongly enhanced couplings to strange 
quarks. Also, since Br(i — > H+b) < 10%, top decays do 
not lead to constraints for our choice of parameters. 

(vii) Additional constraints on the 2HDM can come 
from direct searches for neutral Higgs bosons. The rela- 
tion between the charged and the neutral Higgs bosons, 
however, depends on the details of the Higgs potential. 
For simplicity we assume a Higgs potential such that 
there is one neutral Higgs boson h that has SM-like cou- 
plings to gauge bosons and fermions. Consequently, the 
other two bosons H and A do not couple to gauge bosons 
and have couplings to bottom quarks that are enhanced 
by tan (8/(1 + ebtan/3). While the masses of the Higgs 
bosons are in principle free parameters, we will assume 
Mjj ~ Ma — Mh± to avoid constraints from electroweak 
precision observables. Constraints can arise from neutral 
Higgs boson searches in the H — > tt final state at AT- 
LAS and CMS [77], as well as from searches at D0 and 
CDF in the H -> bb final state [78, 79]. 

In the following, we discuss a benchmark scenario that 
avoids all considered constraints but allows for nonstan- 
dard values for AAqp- We set e„ = i,i t — 0.05i, eg = —I, 
and chose d' d such that i s — 10 -3 . All other are set to 
zero. Fig. 6 shows the regions in the M H ± - tan /? plane 
compatible with the data for AAqp, wa at the la level. 
The green (solid) band is obtained using the expressions 
for the decay amplitude in naive factorization. The blue 
(dashed) band assumes an enhancement of the hadronic 
matrix elements by a factor of 3. The red (dash-dotted) 
region in the upper left corner is excluded by D° - D° 
mixing constraints. 

We observe that sizable NP contributions to the di- 
rect CP asymmetry are only possible for small charged 
Higgs masses and large values of tan/3 ~ 0(100) and 
larger. It is important to note that such large values for 
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FIG. 6. Regions in the M H ± - tan/3 plane compatible with 
the data on Ay4cp, wa at the la level in the discussed 2HDM. 
In a basis with e<j = 0, we set e u = i, it = 0.05i, te = —1 and 
chose t'd such that e s = 10~ 3 . All other a are set to zero. 
The green (solid) band corresponds to the expressions for the 
decay amplitude in naive factorization, blue (dashed) band 
assumes an enhancement of the hadronic matrix elements by 
a factor of 3. The red (dash-dotted) region is excluded by 
D° - D° mixing constraints. 

tan /3 are not in conflict with a requirement of perturba- 
tive Yukawa couplings. Indeed, due to the q, (we find 

0. 03 < ib < 0.09 in the considered region of tan/3) and eg 
factors, the bottom and tau Yukawa remain well below 

1. 4 Extremely large values of tan /3 are not stable under 
radiative corrections. Generically, 1-loop corrections to 
the Higgs potential lead to 1/ tan/3 ~ l/(4-7r) 2 . Values of 
tan /3 ~ few x 100, however, seem a reasonable possibility. 

The shown region of parameter space is not signifi- 
cantly constrained by any of the considered bounds. In 
particular, the bound from Br(Sd — > X s ^) is avoided 
because of the small e t . As stated before, the branch- 
ing ratio of the charged Higgs into rv does not exceed 
10% due to the strongly enhanced couplings to strange 
quarks, and in addition, constraints from top decays are 
satisfied since Br(t -> HH) < 10%. Similarly, we find 
that the branching ratios of the neutral Higgs bosons 
to tt are tiny, O(0.1%), and the stringent bounds from 
ATLAS and CMS [77] are easily avoided. Finally, due to 
the ib factor that controls the size of the bbH and bbA 



4 The phenomenology of generating appropriate bottom and tau 
masses in the very large tan f) regime in the context of supersym- 
mctric models was studied in [80, 81]. We remark however that 
the considered scenario cannot be realized in the MSSM where 
the £i factors are only loop induced and not much larger than 
10~ 2 . 




(a) (b) 

FIG. 7. Example Feynman diagrams that contribute to (a) 
the D — > K + K~ and D — > -k + tt~ decay amplitudes and (b) 
D° - D° mixing in the scalar octet model. 

couplings, the H — > bb searches at Tevatron [78, 79] do 
not lead to constraints. We expect that the collider phe- 
nomenology of this large parameter space will focus on 
a combination of complementary probes, including tests 
for non-SM Yukawa couplings of the light Higgs boson 
and heavy flavor searches for pairs of dijet resonances. 
A complete phenomenological study of this 2HDM with 
MFV is left for future work. 

We stress that the 2HDM with MFV can only signifi- 
cantly affect the direct CP asymmetry in D — >• K + K~ . 
New Physics effects in Aqp(it + tt~) are strongly sup- 
pressed by m d /m s and negligibly small. The considered 
model also does not lead to large nonstandard effects 
in indirect CPV in D° — D° mixing. These statements 
are not necessarily true if the couplings of the Higgses 
radically departure from the MFV ansatz (see [24] for a 
related study). Even though such models will generically 
have flavor changing neutral Higgs interactions, they can 
in principle be made compatible with low energy flavor 
constraints. A detailed study of such setups, however, is 
beyond the scope of this work. 

F. Scalar Octet 

Couplings analogous to the ones in Eq. (32) are also 
possible if the second scalar SU(2)l doublet is a color 
octet instead of a color singlet [82]. We write 

Ant = Y U QUH + Y d QDH^ (44) 
+ X u QT a U<f>% + X d QT a D$ + h.c. 

where H is the SM Higgs and <$jj, a = 1, . . . 8, is the color 
octet. Since a color octet must not aquire a vev, the 
quark masses are entirely provided by the SM Higgs, and 
the Yukawa couplings Y are given by their SM values. 
As in the case of the 2HDM, we will first assume that 
the couplings X have the MFV structure 

X = C Y + C Y Y^Y 

X d = CdY d + C d Y d Y^Y u . (45) 

We can neglect higher powers of the down-type Yukawa 
couplings in the expansion because in the considered sce- 
nario they arc fixed to their small SM values. This auto- 
matically ensures that there are no tree level flavor chang- 
ing interactions of the neutral component of the scalar 
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octet with up quarks. Furthermore, the interactions of 
the charged component of the scalar octet cannot lead to 
tree level contributions to the D — > K + K~ amplitude in 
naive factorization, as the two quarks that couple to the 
color octet cannot hadronize into a kaon or pion, which 
are color singlets. Contributions to the D — > K + K~ 
decay can, however, come from annihilation topologies 
involving the charged octet state. Even though such con- 
tributions are formally suppressed by l/m c , they can be 
sizable in D meson decays. Nonetheless, we now con- 
sider departures from the MFV ansatz for the couplings 
and come back to the MFV framework at the end of the 
section. In particular, we consider the following (small) 
non-MFV entry in the X u coupling matrix in the basis 
where the quark Yukawa couplings are diagonal 



(46) 



with X cu ~ V CS V* S ~ A ~ 0.2. In contrast to the MFV 
interactions, the above term leads to a flavor changing 
neutral current ULCRtfig coupling and, due to SU(2)l in- 
variance, to a 0(1) correction of the d,LCR(f)g coupling 

Ant = CuVcXcu U L T a C R (j>l a 

+ Cuyc(V c * d + X cu V: d ) d L T a c R ^ a + h.c. .(47) 

The tree level exchange of the neutral component of the 
scalar octet, as shown in diagram (a) of Fig. 7, then leads 
to the following contributions to the Wilson coefficients, 



C 



(i) 
si 



l CuCd 



X,, 



4Mi 



(48) 



C 



(i) 

S2 



4 



AC, 



s 

(1) 

T2 



(1) 



Even though the tensor operators do not contribute to 
the D —> K + K~ decay, we include the corresponding 
Wilson coefficients here because they mix with the scalar 
operators under renormalization. 

Tree level contributions to D° — D° mixing from the 
exchange of the neutral complex scalar 4>s can only arise 
if both the flavor changing couplings -O^Cfj^g and clUr^ 
arc present simultaneously. Thus, in our setup, contribu- 
tions to D° — D° mixing first arise at the loop level. The 
dominant contribution comes from a charm quark loop, 
as shown in diagram (b) of Fig. 7, 



C\ 



(2)D 



1 AT? — 2N r + lmt Xl, 



16tt 2 



v 4 8Mi 



(49) 



Similar to the situation discussed in the context of the 
2HDM above, other flavor changing processes like kaon 
mixing or B d — > X s j do not directly constrain or £ d as 
they involve couplings to right-handed top and bottom 
quarks that are proportional to ( t = ( u + (' u yf and 0, = 
Cd + CdVti respectively. Keeping these couplings small to 
avoid constraints from perturbativity and, in particular 
from B d — > V s 7, requires a considerable amount of fine 
tuning. On the other hand, since color octets cannot 



couple to leptons, the constraints from B — > tv and K — > 
[iv are automatically avoided. 

Constraints from electroweak precision observables are 
under control, provided all components of the scalar octet 
have approximately the same mass and are heavier than 
100 GeV [83]. 

There has been extensive work done on collider con- 
straints on scalar octets [84, 85]. Because of the small 
couplings to light quarks in the considered framework, 
single production of scalar octets is strongly suppressed, 
and dijet resonance searches at the hadron colliders give 
no significant constraints in the considered region of pa- 
rameter space. Instead, the main collider constraints 
on the scalar octet doublet come from the dijet pair 
searches at the LHC. We simulate the production of 
scalar octet pairs (including production of neutral pairs 
and charged pairs) using MadGraph 5 and compare the 
obtained cross sections with the bounds set by AT- 
LAS [54] and CMS [55]. We find that the region from 
100 GeV to 200 GeV that is covered by the ATLAS search 
is fully excluded. The CMS search excludes scalar octet 
masses from 320 GeV until approximately 550 GeV for 
13% acceptance (see Subsec. IV C for comments regard- 
ing the 200-320 GeV region). As our scalar octets decay 
predominantly to b quarks, searches for final states with 
multiple b jets at Tevatron can potentially lead to con- 
straints. Using the CDF result [86], the authors of [87] 
find that scalar octet masses of M^ 8 < 200 GeV are ex- 
cluded. Updated results from CDF [79] and D0 [78] 
do not give significantly improved bounds on the corre- 
sponding cross sections. 

In Fig. 8, we show the regions in the M^ s - £„ plane 
that are compatible with the data for A^4cp, wa at the 
lcr level, assuming ( u — ( d /2 real and setting X cu — 
\V cs V* s \e 17 '/ 2 . The green (solid) band corresponds to the 
expressions for the decay amplitude obtained in naive 
factorization, and the blue (dashed) band assumes an en- 
hancement of the hadronic matrix elements by a factor of 
3. The red (dash-dotted) region is excluded by D° — D° 
mixing constraints and the vertical bands are excluded by 
the dijet pair searches at ATLAS (left band) and CMS 
(right band). As expected, constraints from D° — D° 
mixing allow for a sizable A^4cp only for small scalar 
masses. Given the constraints from the dijet pair searches 
at LHC, only a small window with octet masses between 
200 GeV and 320 GeV and large values for C«,d > 20 
remains where significant values for AAcp can be gen- 
erated by the scalar octet with couplings close to the 
MFV ansatz without assuming any enhancement of the 
hadronic matrix elements. Scalar octets heavier than the 
CMS bound of M^ 8 ~ 550 GeV can lead to nonstandard 
values for AAcp only for extremely large Cu.d ^ 60 and 
assuming the hadronic matrix elements to be 3 times the 
naive factorization estimate. The same conclusion holds 
for the pure MFV setup discussed at the beginning of 
the section, provided that the annihilation contribution 
with the charged scalar exchange is of similar size as the 
naive factorization contribution from the flavor changing 
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FIG. 8. Regions in the M^ 8 - <^ u plane compatible with the 
data on AAqp, wa at the lcr level in the scalar octet model, 
assuming ( u = £ d /2 real and setting X uc = | V C sV* s \e 17T ^ 2 . 
The green (solid) band corresponds to the expressions for 
the decay amplitude obtained in nai've factorization, the blue 
(dashed) band assumes an enhancement of the hadronic ma- 
trix elements by a factor of 3. The red (dash-dotted) region is 
excluded by D° - D° mixing constraints. The vertical shaded 
bands are excluded by dijet pair searches at LHC. 



neutral scalar. 

Both in the MFV and the quasi-MFV setup, the scalar 
octet can only significantly affect the direct CP asymme- 
try in D — > K + K~ . New Physics effects in Acp(tt + tt~) 
are strongly suppressed by md/m s and negligibly small. 
As in the case of the 2HDM, this is not necessarily the 
case if the octet couplings are allowed to deviate more 
radically from MFV. We expect collider searches in the 
multi-6 jet final state, as discussed in [84, 85], to provide 
the best sensitivity to this model. 



G. Scalar Diquarks 



Numerous variations of scalar diquarks can have renor- 
malizable couplings to the SM quarks and lead to an in- 
teresting flavor phenomenology [88, 89]. Here, we con- 
sider a diquark that has quantum numbers such that 
it leads to contributions to meson mixing only at the 
loop level, but contributes to the D — > K + K~ and 
D — > tt + tt~ decay amplitudes at tree level in naive factor- 
ization. One such possibility is a scalar that transforms 
as 6 under SU(3), as a singlet under SU(2)l, and has 



hypercharge —1/3. 5 Such a diquark can couple to right- 
handed up and down type quarks, 



Ant — XijUl 



6f 



h.c. 



(50) 



where i,j are flavor indices and the diquark 4>§ is sym- 
metric in the color indices = </>g Q . The couplings 



Xij are a source of flavor violation and in general free 
complex parameters. We remark that the considered di- 
quark can also couple to left-handed quarks. Consider- 
ing couplings to left- and right-handed quarks, however, 
simultaneously induces left-right mixing AF = 2 opera- 
tors, resulting in very severe constraints from D and K 
meson mixing. Therefore, we concentrate on couplings 
to the right-handed quarks only. 

The tree level contributions of the diquark to the D — > 
K + K~ and D — > tt + tt~ decay amplitudes are shown in 
diagram (a) of Fig. 9 and read 



(i)p 



Y Y* 

2 " TgmT 



(51) 



Since QCD is parity conserving, the matrix element of 
the operators &i are the same as the SM operators 
0^ p ', and no strong phase difference occurs between 
them. Correspondingly, c[^ p does not contribute to the 
direct CP asymmetries (see Eq. (16)). Due to the dif- 
ferent color structure of the operators 0% , however, 
a different strong phase can be expected in their ma- 
trix elements and a non-zero contribution to the direct 
CP asymmetries can be generated by the weak phase in 



a 



oc X C pX up . 



Contributions to D° — D° mixing are first induced at 
the loop level (see diagram (b) in Fig. 9) 



(2)D 



1 N c + 3 
16tt 2 32M 2 



^ X u P x c 

\p—d,s 



(52) 



where, for simplicity, we restrict ourselves to couplings 
of the diquark to the first two generations of quarks and 
expand the result in small quark masses. 

We note that even with small couplings to quarks 
< 0.1, sizable effects in flavor observables can be gen- 
erated. Consequently, dijet searches at hadron colliders 
do not lead to relevant constraints. Important collider 
constraints on diquarks, however, come from the dijet 
pair searches at LHC. Using MadGraph 5 to simulate di- 
quark pair production, we find that the ATLAS search 
excludes the region from 100 GeV to 200 GeV. The CMS 
search excludes diquarks with masses from 320 GeV until 
1000 GeV assuming an acceptance of 13%. 



5 All other possibilities are scalars that can also have lepton num- 
ber violating couplings. One example is a SU(3) triplet, 5C/(2)^ 
singlet with hypercharge —1/3 which corresponds to a right- 
handed down squark with i?-parity violating couplings. 
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FIG. 9. Example Feynman diagrams that contribute to (a) the D — » K + K and D — > tt + it decay amplitudes, (b) D° — D° 
mixing and (c) K — K mixing in the discussed diquark model. 
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FIG. 10. Regions in the - \X cs X^ g \ plane compatible 

with the data on A^cp, wa at the la level in the diquark 
model, setting Arg(X cs X^ s ) = 7r/2. The green (solid) band 
corresponds to the expressions for the decay amplitude ob- 
tained in naive factorization, the blue (dashed) band assumes 
an enhancement of the hadronic matrix elements by a factor 
of 3. The red (dash-dotted) region is excluded by D° — D° 
constraints. The vertical shaded bands are excluded by dijet 
pair searches at LHC. 



In Fig. 10, we show regions in the M^ 6 - \X CS X* S | plane 
compatible with the data for AAgp, wa at the la level, 
considering only X cs and X us to be non-zero and setting 
Arg(X cs X* s ) = 7r/2. The green (solid) band corresponds 
to the expression for the decay amplitude in naive fac- 
torization. The blue (dashed) band assumes an enhance- 
ment of the hadronic matrix elements by a factor of 3. 
The red (dash-dotted) region is excluded by the D° — D° 
mixing constraints. The vertical bands are excluded by 
the dijet pair searches at ATLAS (left band) and CMS 
(right band). Since only strange quark couplings are con- 
sidered, AAqp is entirely generated by Acp(K + K~) in 
this setup. Despite the fact that the D — > K + K~ decay 
amplitude arises already at the tree level while D° — D° 
mixing is only induced at the loop level, constraints from 



D° — D° mixing combined with the dijet pair searches 
exclude a sizable AAqp unless the hadronic matrix el- 
ements are enhanced by a factor more than 3. This is 
due to the fact that the contributions to AAqp coming 

~ ( l)s 

from the operator 0\ are color suppressed and neither 
enhanced by RG effects nor by chiral factors. Consider- 
ing only the couplings X cc i and X u( i leads to an analo- 
gous situation, but in this case AAqp entirely stems from 
A C p(7r + 7r~). 

If all 4 couplings X cs , X us , X cd and X ud are present, 
the D° — D° mixing constraint can be strongly reduced by 
assuming a GIM-like mechanism, i.e. assuming the cou- 
pling matrix X to be (approximately) unitary. This hap- 
pens, for example, if X is (approximately) proportional 
to unity in the weak eigenstate basis for the quarks. In 
that case, however, Acp(tt + tt~) ~ Aqp{K + K~) and the 
difference between the two, AAqp, is very small. 

In addition, also loop contributions to kaon mixing are 
induced if all four couplings are present. Apart from a 
contribution analogous to Eq. (52), mixed <pg - W loops 
(shown in diagram (c) of Fig. 9) can become very impor- 
tant. For the latter, we find 

cf^ K = cP K = y|^- (v cs v* d ) (X cs X* d ) 

1 m 2 c ( rn 2 c \ 

where we expand the result in the charm quark mass 
and keep only the leading term that is enhanced by a 
large logarithm. Despite the suppression by m^/M^, 
these contributions are very relevant because the matrix 
elements of 0^ K and 0^ K are strongly chirally en- 
hanced and enhanced by renormalization group effects. 6 
Furthermore, the simultaneous presence of X us and X ud 
also leads to tree level contributions to e'/e. 

Even after tuning phases in order to minimize con- 
straints from CPV in kaon mixing and e'/e, we do not 



6 The analogous contribution to D° — D° mixing is suppressed by 
the strange quark mass and has a much smaller chiral enhance- 
ment from the matrix elements: it is therefore negligible. 
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find any region in parameter space of the considered di- 
quark model in which a sizable AAqp is compatible with 
all constraints. 



V. NEW PHYSICS CONTRIBUTIONS TO 
GLUON PENGUINS 

We now consider possibilities for new physics in the 
gluon penguin diagrams at the 1-loop level. We con- 
sider models with new fermion and scalar fields, where 
the fermion is Dirac or Majorana and with or without 
GIM suppression in Subsec. VA and Subsec. VB, re- 
spectively. A scenario with chirally enhanced magnetic 
penguins is considered in Subsec. V C. As discussed in the 
Introduction, the chiral enhancement of the decay ampli- 
tudes allows such a scenario to be the least constrained 
by D° — D° mixing. 



Fermion + Scalar Loop without GIM 
Mechanism 



We add one heavy fermion \ an d one heavy scalar <fi 
to the SM that couple with right-handed up type quarks 



X u urx4> + X c c R x<fi + Xtt R x(l> 



h.c. 



(54) 



The couplings X± necessarily violate flavor and are free 
complex parameters. In the following we consider the 
case where the new particles are SU{2)l singlets. We as- 
sume the fermion is electrically neutral and a SU(3) sin- 
glet and the scalar is a SU (3) triplet. Different quantum 
number assignments do not lead to qualitatively differ- 
ent results concerning the D — > K + K~ and D — > 7r + 7r~ 
decay amplitudes. We do not consider new fields with 
flavor breaking couplings to left handed quarks to avoid 
constraints from kaon mixing and e'/e. 

The Lagrangian considered in this framework incorpo- 
rates an accidental Z2 symmetry where \ an( l 4> arc °dd 
under the symmetry while the SM fields are even. In 
the considered minimal framework, the lightest of x an d 
<j> is therefore stable and becomes the dark matter. If 
the colored scalar is lighter than the fermion, additional 
interactions have to be introduced to prevent the scalar 
from being absolutely stable. If these additional inter- 
actions are weak enough such that the scalar remains 
stable on detector scales, bounds from searches for long- 
lived particles apply [90, 91]. If instead the scalar decays 
promptly, then it is highly model dependent whether or 
not the considered scenario is excluded by current collider 
searches. 

In the case where the neutral fermion is the lighter 
particle, SUSY searches with jets + $ T at the Teva- 
tron [92, 93] and LHC [94, 95] provide constraints to the 
model. In particular, both ATLAS [94] and CMS [95] 
give bounds on the production cross section for a simpli- 
fied model that contains the first and second squark gen- 
erations and a neutralino lightest supersymmetry parti- 



cle (LSP). We simulate the production cross section using 
Prospino 2 . 1 [96] and set all the superpartner masses 
to be 4.5 TeV except the relevant squarks and the LSP, 
following [94]. The obtained cross section is then scaled 
by 1/8 since in our minimal scenario there is only one 
colored scalar instead of eight. We find that neither the 
ATLAS nor the CMS searches currently put bounds on 
our model. 

The discussion of the collider constraints is the same 
regardless of whether the neutral fermion is a Majorana 
or a Dirac particle. On the other hand, the nature of 
the fermion plays a crucial role in the low energy phe- 
nomenology of the considered model. 



1. Dirac Fermion 

The interactions specified in Eq. (54) lead to 1-loop 
gluon penguin contributions to the Wilson coefficients 
of the AF = 1 effective Hamiltonian that are shown in 
diagram (a) of Fig. 11 and lead to 

CP = ^X u X* c7 ^p(z) , 



An 



8Mi 



Mi) _ _ _ J_r-(i) _ _ J_r<(! 



Nr. 



a 



(1) _ 



X U X* n g(z) 

c 4M? 



(55) 



The loop functions depend on the ratio of the fermion 
and scalar masses z = M^/M^. We find p(l) = -1/24, 
g(l) = 1/48. Their full analytical expressions can be 
found in Appendix C. The relation between the Wilson 

coefficients 5 6 of the QCD penguin operators is uni- 
versal for all models where they are induced by gluon 
penguins. 

The couplings in Eq. (54) also lead, in principle, to 
1-loop box contributions to the AF = 1 effective Hamil- 
tonian. These contributions, however, are strongly sup- 
pressed by an additional factor of \X U \ 2 and are therefore 
negligible. 

Finally, 1-loop box contributions to D° — D° mixing 
are also induced (see diagram (b) of Fig. 11) 



-^n = {X u X* c f 1 



16tt 2 Ml 8 



/(*), 



(56) 



with /(l) = —1/3. The analytical expression for / is 
given in Appendix C. 

We find that D° — D° mixing constraints exclude NP 
contributions to Aj4cp as l a rge as the world average even 
if we allow for an enhancement factor of 3 in the hadronic 
matrix elements. 



2. Majorana Fermion 

We find qualitatively different results for the Majorana 
fermion. While the gluon penguin contributions to the 
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FIG. 11. Example Feynman diagrams that contribute to (a) the D — s> K + K and D — > tt + tt decay amplitudes and (b and c) 
D° - D° mixing in the discussed models with fermion + scalar loops. 
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FIG. 12. The M$ - M x plane in the models where a gluon penguin is induced by a Majorana fermion - scalar loop without GIM 
mechanism (left) and with GIM mechanism (right). The blue contours indicate the (left) = 0.06, ...0.14 or (right) 

\5 CU \ ~ 0.05, . . . 0.25 values in agreement with the measured AAcp, wa, assuming maximal phases of n/2. An enhancement 
of the hadronic matrix elements by a factor of 3 is assumed in both plots. The red (dash-dotted) region is excluded by the 
D° - D° mixing constraints. In the region above the black (solid diagonal) line, the colored scalar is lighter than the fermion. 



decay amplitudes are identical for Majorana and Dirac 
fermions, the box contributions to D° — D° mixing dif- 
fer. In the case of a Majorana fermion, "crossed" box 
diagrams also exist (see diagram (c) of Fig. 11). One 
finds 

6r=^^(i/w+i/w), < 5 7> 

with /(l) = 1/6. The analytical expression for / can be 
found in Appendix C. 

The left plot in Fig. 12 shows the - M x plane in 
the considered scenario. The blue contours labeled with 
0.06, ... 0.14 show the \X U X* | values that are required to 
be in agreement with the measured AAqp, wa at the la 
level, assuming Axg(X u X*) = n/2. An enhancement of 
the hadronic matrix elements by a factor of 3 is assumed. 
The red (dash-dotted) region is excluded by the D° — D° 



mixing constraints. In the region above the black (solid 
diagonal) line shaded in dark gray, the colored scalar is 
lighter than the fermion and would be a stable colored 
particle: for this region of parameter space, the model 
must be extended to be viable. 

We observe that if the Majorana and scalar masses are 
equal, Mtp — M x , the constraint from D°~D° mixing can 
be avoided. Indeed, for such a ratio of masses, the box 
and crossed box contributions to the mixing amplitude 
cancel. The exact mass ratio where such a cancellation 
occurs depends on the quantum numbers of the scalar 
and the Majorana fermion. For example, in the well know 
SUSY case of squark - gluino boxes, the cancellation oc- 
curs for Mg ~ 1.6 x Mq (see [97] for a recent discussion.). 
Because of this cancellation, there exists a narrow region 
of parameter space where the NP contribution to AA<jp 
can explain the measured value. The same is in principle 
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true without assuming any enhancement of the hadronic 
matrix elements. In such a case, however the D° — D° 
mixing constraints hardly allow any deviation from the 



M = M x line. 



B. Fermion + Scalar Loop with GIM Mechanism 

We consider now a framework that allows the incor- 
poration of a GIM-like mechanism. In addition to the 
heavy fermion, we introduce 3 scalars, <fi u ,c,t, which are 
scalar partners of the 3 right-handed up-type quarks. If 
the scalars couple universally to quarks and are approxi- 
mately degenerate in mass, summing over the three scalar 
loops strongly suppresses flavor changing neutral current 
processes. As discussed in [98] for this setup, AF = 2 
processes are more strongly GIM suppressed with respect 
to AF — 1 processes and we expect more room for NP 
in the decay amplitudes. 

We consider the following flavor universal couplings 
between the quarks and scalars 

Ant = 94>u R X<t>u + g^CRX^c + 9^i R x<t>t + h.c. . (58) 

In the basis where the quark masses are diagonal, we 
write the mass matrix for the scalars as 



Ml = Mfi 



+ Ml8 



(59) 



where Sij <C 1 are a source of flavor violation. Such a 
setup is, for example, realized in the MSSM by gluinos 
and right-handed up squarks. Here, we focus on the 
case where the new particles are SU(2)l singlets, the 
scalars are SU(3) triplets, and the fermion is a SU(3) 
singlet. Different quantum number assignments do not 
lead to qualitatively different results concerning the D — > 
K + K~ and D —> it + tt~ decay amplitudes. 

Similar to the scenario without the GIM mechanism 
discussed above, the framework considered here is only 
mildly constrained by collider searches. The difference 
with respect to the case without GIM mechanism is 
that now there are three scalars, and the production 
cross section for scalar pair production is therefore three 
times larger than in Subsec. VA. On the other hand, 
one of the scalars decays predominantly into top quarks 
and the constraint from squark pair production with di- 
rect decays into light quarks and the lightest neutralino 
from [94, 95] cannot be applied for this scalar. Analo- 
gously to the previous scenario without GIM, we scale 
the production cross section obtained by Prospino 2.1 
by the ratio of the number of scalars in the model, i.e. 
a 2/8 scaling. We hnd that the ATLAS search [94] does 
not put any bound on our model. The CMS search [95] 
on the other hand excludes a small corner of parameter 
space with scalar masses between 300 GeV and 350 GeV 
and fermion masses below ~ 100 GeV. This holds regard- 
less of whether the fermion is Majorana or Dirac. The 
flavor phenomenology is, however, qualitatively different 
for the two cases. 



1. Dirac Fermion 

The 1-loop gluon penguin contributions to the decay 
amplitudes are 



a 



(i) 



4tt uc 8M2 



=,(1) 



1 X. 



(1) 



1 = 



(1) 



Nr. 



Nr. 



9% 



(60) 



where S uc = S* u is a complex dimensionless parameter 
of flavor violation as defined in Eq. (59). For the loop 
functions, we find P(l) = 1/30, G(l) = -1/80. Their 
analytical expressions can be found in the appendix. 

Now, 1-loop box contributions to the AF = 1 effec- 
tive Hamiltonian are not suppressed by additional small 
mixing angles. For the minimal set of couplings defined 
in Eq. (58), we find 



1 X 



(1) 



16tt 2 YIMl 



B(z) 



(61) 



where B(l) = 1/48. Its analytical expression can be 
found in Appendix C. We remark that in naive factor- 
ization, this combination of Wilson coefficients does not 
contribute directly to the D — > K + K~ and D — > ir + ir~ 
decays that have either strange or down quarks in the 
final state. 7 Through renormalization group running, 
however, the other QCD penguin operators (in particu- 
lar Oq ) are induced and lead to non-zero contributions. 
Although such contributions are subleading, they are in- 
cluded in our numerical analysis. 

Finally, the 1-loop box contributions to D°—D° mixing 
are given by 



(2)D 



16tt 2 M 



9 \\m 



(62) 



with F(l) = —1/30. The analytical expression for F can 
be found in Appendix C. 

Similar to the situation without a GIM mechanism, 
we find that the gluon penguins that are induced by a 
Dirac fermion loop cannot viably account for the mea- 
sured AAcp t wAi even if an enhancement of the hadronic 
matrix elements by a factor of 3 is assumed. 



2. Majorana Fermion 

The gluon penguin contributions to the decay ampli- 
tudes, that are induced by a Majorana fermion - scalar 



7 They do contribute in QCD factorization via annihilation dia- 
grams that we do not consider here. 
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loop, are identical to the Dirac case, but the box contri- 
butions to the decay amplitudes and to D° — D a mixing 
differ. Adding the "crossed" boxes, we find 
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(i) 



1 /=V(l) 

2° 9 



(b(z)--B(z) ) . ((«) 
16tt 2 12M? V 2 v ' 1 



(2)D 



si 

16tt 2 Ml 



9% 



(64) 



with B(l) = -1/12, F(l) = 1/20. Their analytical ex- 
pressions can be found in Appendix C. 

The corresponding situation in the M^ - M x plane is 
shown in the right plot in Fig. 12, assuming an enhance- 
ment of the hadronic matrix elements by a factor of 3. 
The blue contours labeled with 0.05, . . . 0.25 show the 
\5 UC \ that are required for agreement with the measured 
AAcp, wa at the la level, setting Arg((5 uc ) = ir/2. Now, 
the ratio of the scalar to Majorana mass where the box 
contributions to meson mixing cancel is M^ ~ 2.3 x M x . 
We observe that constraints from D° — D° mixing leave 
more regions of parameter space open, as expected. Still, 
the measured value for A^4cp can be explained by NP 
contributions only in a narrow region along the M^ ~ 
2.3 x M x line. 



C. Chirally Enhanced Magnetic Penguins 

Finally, we discuss a setup that leads to chirally en- 
hanced chromomagnetic penguins. We introduce a Majo- 
rana fermion \ that is singlet under the SM gauge group, 
as well as scalar partners, <j) q and <p u , to the left-handed 
quark doublets Q and the right-handed up-type quarks 
U, respectively. We consider the following interactions 
among those degrees of freedom 



Ant, mass — X u Ux<l}u + X q Q X (f> q + A<j)* u 4> q H + h.c. 



+ml\<j) u \ 2 + m 2 q \0 q \ 2 + M xX X 



(65) 



We set the coupling matrices X u and X q to be universal 
in flavor space: X q — g q I, X u = g u l. We further assume 
the mass matrices for the scalars m 2 and m 2 are universal 
and, for simplicity, also equal: m 2 — m 2 = Mjl. The 
only new source of flavor violation is then the trilinear 
coupling A. After electroweak symmetry breaking, the 
trilinear coupling A leads to mixing between the 4> u and 
the isospin +1/2 component of <p q . This setup resembles 
to a large extent the MSSM with flavor changing trilinear 
couplings in the up-squark sector as discussed in [18, 26]. 

Both the elements A cu and A uc can lead to c — > u 
transitions. Expanding the contribution to the Wilson 
coefficients in A cu v/M 2 < 1, we find 



a 



a 



(i) 

(1) 
83 



Mi 



AMI 1 ' 



g q 9u v ^* UC ^ X -G(z) 



M 2 



AMI 



(66) 



where z = M 2 /M 2 and (5(1) = -1/24. The full analyti- 
cal expression for the loop function G can be found in Ap- 
pendix C. We highlight the expected chiral enhancement 
of the magnetic penguins by a factor of v/m c . 

Box diagrams with the new fermion and the scalars 
also lead to contributions to D°—D° mixing. Considering 
only the A cu coupling, we find 



C (2)D 



2 2 

9 q 9u 



1 (A cuV ) 2 l p{zh 



16tt 2 Ml M 
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where the loop function F is the same as in Subsec. VB. 
The A uc coupling leads to the analogous contribution 
to the coefficient C^ D . If both couplings are present 



simultaneously, contributions to C^ D and C y ^ ,1J are also 
generated. 

In Fig. 13, we show regions in the M^ - g u g q \A cu \/M^ 
plane compatible with the data on A^4cp. wa at the la 
level, considering A cu as the only source of flavor viola- 
tion and setting Arg(A cu ) = tt/3. In the plot on the left 
(right) we set M x = 100 (250) GeV. The green (solid) 
band corresponds to the expression for the decay am- 
plitude in naive factorization. The blue (dashed) band 
assumes an enhancement of the hadronic matrix elements 
by a factor of 3. The red (dash dotted) region is excluded 
by the D° — D mixing constraints. The gray shaded re- 
gions are excluded by SUSY searches with jets + $t- 
Results from ATLAS [94] and CMS [95] indicate that 
scalar masses up to 675 GeV are excluded for a fermion 
mass of 100 GeV (see left plot). On the other hand, for 
a fermion mass of 250 GeV, only a small region around 
500 GeV is excluded by CMS data alone. For heavier 
fermion masses, the full range of scalar masses is allowed 
by present collider constraints. 

Because of the chiral enhancement factor v / m c in the 
decay amplitudes, D° — D° mixing constraints allow for 
a large AAqp in the considered setup. Even for scalar 
masses of 1 TeV and larger, NP contributions to D° — D° 
mixing are more than an order of magnitude below the 
experimental constraints for NP contributions to AAcp 
that agree with the experimental value. 



,(2)13 



VI. CONCLUSIONS 

The LHCb measurement of the difference in the time 
dependent CP asymmetries in the singly Cabibbo sup- 
pressed D — > K + K~ and D — > ir + ir~ decays, AAqp, is 
the first evidence for charm CP violation. Although there 
are large uncertainties in the SM prediction, the measure- 
ment could indicate New Physics, and NP interpretations 
are nevertheless motivated and exciting. In this paper, 
we studied the effect of NP degrees of freedom for non- 
standard direct CPV in the D -> K + K~ and D — > ir + ir~ 
decays while also considering constraints both from low 
and high energy experiments. 

As is shown in [21], models that give rise to chirally 
enhanced chromomagnetic c — > u penguin operators are 
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FIG. 13. Regions in the M$ - g u gq\A cu \/ plane compatible with the data for Aj4cp, wa at the la level in the model with 
chirally enhanced gluon penguins, setting Kxg(A cu ) = 7r/3. In the left (right) plot M x — 100 (250) GeV. The green (solid) band 
corresponds to the expression for the decay amplitude in naive factorization, the blue (dashed) band assumes an enhancement 
by a factor of 3. The red (dash-dotted) region in the upper left corners is excluded by the D° — D° mixing constraint. The 
dark gray vertical region is excluded by jets + $ T searches at LHC. 



the least constrained by low energy data and can easily 
accommodate the large AAcp value measured by LHCb. 
The most prominent examples for such models are super- 
symmetric scenarios as discussed in [18] and very recently 
in [26]. We studied the chirally enhanced chromomag- 
netic penguins in the framework of a simplified model 
that contains scalar partners of the left- and right-handed 
up-type quarks as well as a Majorana fermion. We con- 
firm that low energy observables, in particular D° — D° 
mixing, as well as collider searches do not significantly 
constrain the model's parameter space that leads to a 
sizable A^4cp- 

Models that contribute to the D meson decays through 
four fermion operators are generically expected to be 
strongly constrained by D° — D° mixing data [21]. Non- 
standard effects in the decays are only possible if the 
new degrees of freedom mediating the c — > u transition 
are very light. In this work, we quantified this statement 
through a systematic study of models with a minimal set 
of new degrees of freedom giving rise to four fermion op- 
erators both at the tree and the loop level. In summary, 
we find: 

• Flavor changing couplings of the the SM Z boson 
can induce a AAcp as large as the observed value 
if the NP phase is moderately tuned to avoid con- 
straints from indirect CP violation in D° — D a mix- 
ing. 

• A Z' that mediates the c —> u transition at tree 
level cannot account for the observed A^4cp due to 



the combined constraints from D° — D° mixing and 
dijet searches. 

• A heavy gluon with a flavor changing tree level 
c — » u coupling and with a mass of 200 GeV 
< Mg> ;$ 320 GeV cannot fully be excluded as 
NP explanation for the measured AAqp if a mod- 
erate enhancement of the hadronic matrix elements 
is allowed. 

• In a 2HDM with MFV, there exist regions of pa- 
rameter space that can lead to a sizable A^cp. 
They are characterized by light charged Higgs 
masses and strongly enhanced couplings of the right 
handed strange quark to the charged Higgs, with 
respect to its SM Yukawa coupling. Avoiding con- 
straints from perturbativity and B — > X s j, how- 
ever, requires a considerable amount of fine tuning. 

• Scalar octets can also induce large nonstandard ef- 
fects in a mass window 200 GeV < M 08 < 320 GeV 
that is left open by current collider searches. The 
viable parameter space is analogous to the 2HDM 
with MFV model and appears to be, to some ex- 
tent, fine tuned. 

• The scalar diquark model we consider is ruled out 
by D° — D° mixing as a NP explanation of A Acp ■ 

• The minimal models that induce NP effects in the 
D meson decays through loops with Dirac fermions 
and scalars are strongly constrained by D° — D° 
mixing data and cannot give rise to a sizable A Acp ■ 
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• If Majorana fermions and scalars appear in the 
loops, then D° — D° mixing constraints can be 
avoided for a particular ratio of Majorana and 
scalar masses that depends on the exact quantum 
number assignment for the particles. Correspond- 
ingly, in such models there exist regions of param- 
eter space that lead to large A^4cp m agreement 
with the data. 

We note that our results are robust, since changes to 
the central value of AAqp, new direct search constraints, 
and the enhancement or suppression of flavor bounds 
coming from additional field content can be readily ap- 
plied to our minimal models and our derived Wilson co- 
efficients. 

As we showed, the New Physics parameter space fa- 
vored for an explanation of the LHCb evidence for charm 
CP violation is largely within the current reach of various 
direct searches at the LHC. Thus the well-known comple- 
mentarity between low energy flavor measurements and 
high energy direct probes may prove fruitful as we con- 
tinue to search for New Physics. Our work emphasizes 
this synergy by presenting a broad study of minimal New 
Physics models, discussing both their effects on low en- 
ergy flavor observables as well as their high energy col- 
lider signatures. We have demonstrated that a number 
of intriguing New Physics models can viably explain the 
large AAcp measurement, and we have concretely iso- 
lated the interesting parameter spaces of such models 
which must now be searched directly. 
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Appendix A: Hadronic Matrix Elements in Naive 
Factorization 

To evaluate the hadronic matrix elements of the oper- 
ators in the AF = 1 effective Hamiltonian in Eq. (14), 



we use naive factorization 



(K+K-\(ur lS )(sT 2 c)\D°) 
^(K+\(uT lS )\0)(K-\(sT 2 c)\D°) , (Al) 

where Ti represent the various Dirac and color structures. 
In this approximation, it is straightforward to evaluate 
the hadronic matrix elements (Oi) = (K + K~\Oi\D°) 

(0«) = N c {0i\= N c {0^) = (OP) (A2) 
= -27V C (0«) = -2<0#) , 

J-( W) = N c {OP) = (0™) = -27V C (0«) 
Xf 

= -2(0«) = -8<0&) = -8iV c (0«) , 
(0$) = (0$) = , 

where N c = 3 is the number of colors and X/ is the appro- 
priate chiral factor from Eq. (17). Using QCD factoriza- 
tion methods, the matrix element of the chromomagnetic 
operator is [18] 

= -S^P + X)(0<"). (A3) 

As QCD conserves parity, the matrix elements of the 
chirality flipped operators are identical to the ones 
shown above. 



Appendix B: Anomalous Dimensions 

For completeness we collect here all the anomalous di- 
mensions of the AF = 1 operators that are required for 
our analysis. The LO anomalous dimension matrix that 
governs the running and mixing of the current-current 
operators ', the QCD penguin operators 6 and 
the QED penguin operators 10 is given by (see 

e-g. [99]) 
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where N c — 3 is the number of colors, / is the number 
of active quark flavors, and u and d are the numbers of 
active up- and down-type quarks, respectively. 

Leading order running of the chromomagnetic operator 

°8ff is s iven b y 



mi - 8 



(B2) 



For the LO anomalous dimension matrix responsible for 
the running and mixing of the scalar and tensor operators 
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which agrees with [31] once the different conventions for 
the operators and the <7 M „ matrix are taken into account. 
At leading order and in the limit of massless down and 
strange quarks, the scalar and tensor operators do not 
mix into other operators. 

As QCD conserves parity, the anomalous dimensions 
for the chirality flipped operators O 
the ones shown above. 
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are identical to 



Appendix C: Loop Functions 

The loop functions hi and /12 appear in the charged 
Higgs contributions to kaon mixing in the 2HDM in Sub- 
sec. IV E 
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The loop functions p, g, f and / appear in the discus- 
sion of the model with fermion - scalar loops without a 
GIM mechanism in Subsec. VA. The functions p and g 
occur in the expressions for the gluon penguin contribu- 
tions to the D meson decay amplitudes. The function / 
comes from the evaluation of a box diagram contributing 
to D° — D° mixing and / comes from the corresponding 
crossed box diagram 
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The loop functions P, G, F and F are the analogues 
to p, g, f and / in the fermion - scalar model with a 
GIM mechanism discussed in Subsec. VB. The functions 
B and B appear in the box and crossed box contributions 
to the D meson decay amplitudes of that framework: 
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Finally, the loop function G appears in the expression 
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for the chromomagnetic penguin loop in Subsec. V C 

n, \ 1 + 5^ z(2 + z) 
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